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ABSTRACT 

Use of a low level pilot tone has been shown to eliminate the 
error floor in fading channels. This paper demonstrates that 
non- idealities in the receiver's pilot tone filter cause reappearance 
of the error floor. It also presents the BER in closed form, in 
contrast to the multidimensional numerical integration of previous 
work. 

INTRODUCTION 

Recently, a number of papers [e.g. Davarian, 1986, 1987; Simon 
1986; LaRosa & Citron, 1987] have resurrected the use of a pilot tone 
to allow coherent reception in fading channels . The receiver extracts 
the pilot with a pilot bandpass filter (PBPF) , and uses the result as 
a phase reference. The obvious costs are the non- constant envelope, 
the additional complexity at the receiver, and the fact that the pilot 
robs power from the data carrying signal. Thus we expect to see 
little or no improvement in the noise limited region; however, when 
the fade rate is high, TCT proves itself, by lowering or eliminating 
the error floor typical of other modulation methods. 

Analyses to date have implicitly or explicitly assumed that TCT 
eliminates the error floor. We show here, though, that nonidealities 
In the PBPF lead, not just to a degraded BER, but to a nonzero error 
floor. Moreover, previous studies have calculated the BER by 
numerical integration, often in multiple dimensions. We present 
here a closed form solution which includes the effect of the PBPF 
nonidealities and a frequency offset error. 

In the interest of brevity, this paper deals only with the issue 
of irreducible error rate, because it has not been addressed by other 
authors. However, the method extends in a straighforward way to 
include the effect of noise, and a full analysis will be published 
elsewhere . 

SYSTEM MODEL 

General Characteristics 

Figure 1 is the complex envelope representation of the 
transmitted and received signals. The transmitted power is split 
between the data bearing signal and a low level pilot tone in such a 
manner that they do not interfere with each other. For simplicity, we 
will assume Manchester coding with the pilot in the center channel 
null, although other methods are possible; for example, the pilot 
could be located beside the data signal. The analysis below is 
restricted to two-level eyes; for convenience we assume BPSK, though 
the results are also valid for QPSK or OQPSK. 
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Data is detected with a filter matched to the transmitted pulse 
p(t). Samples at the filter output are denoted by u(k) . 

The PBPF , represented by the low pass equivalent filter H^(f) , 

has unity dc gain; ideally, it is a rectangular lowpass filter with 
bandwidth B /2, made equal to the sum of frequency offset and Doppler. 

Since the PBPF has a greater delay than the wider bandwidth matched 
filter, we assume that there is delay equalization in the detector 
branch! The PBPF output is v(k) . Feedforward phase correction is 
accomplished simply by making decisions on the basis of the 
phase -corrected decision variable: 

d(kT) - Re[u(kT) v*(kT)] (1) 

The shape and bandwidth of the PBPF have a major effect on 
performance. If it is too narrow, then frequency offset and Doppler 
spread will prevent the filter output from following the channel 
fluctuations, leading to degraded performance and a nonzero error 
floor. If it is too wide, it will admit too much channel noise and 
self noise from the data signal, which again leads to degraded 
performance. A related issue is the fact that any phase distortion in 
the PBPF also results in decorrelation and a non zero error floor. 

We will assume that there is no additive noise in the channel; as 
noted earlier, our main interest here is the irreducible BER. 


The Fading Channel 

We assume Rice fading, in which the ratio of specular to diffuse 
power is K. These components of the complex gain, when each 


normalized to unit power, have power spectral densities S (f) and 

S (f) , respectively. Although our general solution does not rest on 
gd 1 

a particular spectrum for the diffuse component, the examples will use 
a spectrum and autocorrelation function characteristic of isotropic 
scattering: 


V f) - ; < f D ' f2) ' 1/2 ; V° ‘ V 2 'V> <2) 

where J () is the Bessel function of first kind and order zero and f^ 
o 

is the Doppler frequency. For the specular component, we have: 

S <f) = fi(f-f L ) (3) 

where is its Doppler shift. Note that | f^ | — fjy - We also all° w a 
frequency offset error f between transmit and receive oscillators . 


Signal Moments 

The signal moments and their ratios are required for calculation 
of BER. Since we have assumed no noise, the variance in u(k) and v(k) 
is due to the diffuse component, and the mean is determined by the 
specular component. We can show that the samples at the output of the 
matched filter have signal to variance ratio: 
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Similarly, at the pilot tone filter output: 



> 
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(5) 


l"v' 2/ % - 2K l H p< £ o +f L , | Z ( J S gd <f ' f o )|H p <£)|2df ) 

Next, we calculate the normalized inner product of the means of u(k) 
and v(k) : 

* * f r~ 9 l -1/2 

VA*» " 2K H p (f o +f L ) [ J S gd <f - f o ) I H p<« I « J <6> 

Finally, we calculate the correlation coefficient of u(k) and v(k) : 

JV f - f o )H p< f>df (l S gd <f ‘ f o )|H p <f)|2df ) 1/2 <7> 

DETECTION ERROR PROBABILITY 
General Solution 

The basic observation is that u(k) and v(k) are jointly Gaussian, non 
zero mean, complex random variables. The decision variable d(k) is 
then of a form familiar from studies of incoherent detection; we make 
an error if it is negative, assuming b(k) is positive. A concise 
derivation of the probability of this event can be found in Schwartz 
et al, 1966, sec8.2. Define the quantitities a and b : 



and the quantity C: 

f 1,2 -w 1/2 

C “ [C 1 * + p x\ / 2 C X ‘ p i) (9) 


where p - + p Then the probability of bit error is: 

P e “ Q(a,b) - C I Q (ab) exp( - (a 2 +b 2 )/2) (10) 

where Q(a,b) is the Marcum Q-functi®n, and I q () is the modified Bessel 

function of first kind and order 0. Techniques for evaluating the Q 
function can be found in Bird & George, 1981. 

In the case of Rayleigh fading, K - 0, the parameters a; and b are 
zero, and the BER simplifies to: 

f o 1/2 1 1 > 2 

P eR - 1 - C - [t 1 - 'i) - /- r J / 2(1 - p[) .... (U) 

Moreover, if the frequency response H^(f) is real, then = 0, and: 

P eR = (l-p)/2 (12) 

Equations (11) and (12) illustrate the requirement for extremely high 
correlation between the outputs of the matched filter and the PBPF. 

We see that filters other than the ideal rectangular lowpass produce a 
correlation coefficient less than 1 in magnitude, and therefore an 
irreducible error rate. 

APPLICATIONS 

Delay Mismatch 
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Earlier we noted that the greater delay of the narrow PBPF meant 
that delay equalization is required in the matched filter branch. In 
an implementation based on DSP, this is unlikely to be a problem. In 
an analog implementation, though, it is possible to have a residual 
delay mismatch. We can incorporate this effect by representing an 
otherwise perfect PBPF as: 

H (f) = exp( - j 2jrfr ) (I 3 ) 

P 

where the delay mismatch r can be positive or negative. We now have: 

p - R At) exp(-j2wf t) ( 14 ) 

gd 0 



K 


1 + RgSin(2ir(f o +f^)r) T cos(2ir(f o +f^)r) Vl - Rg 

(1 - 


C- 0.5 [ Vl ■ R* + R c ] (l - Rg r 1/2 <15) 

where for notational convenience we have defined the quantities . 

R C " R gd (r) cos ( 2wf o T) : R s “ R gd (r) sin ( 2jrf 0 r) (16) 

Once again, the case of Rayleigh fading is particularly simple. 

We set K - 0 and find, as in (11), that the error floor is: 

P eR- 0.5 .R c ) Cl-R 2 r V2 <”> 

Figure 2 illustrates the effect of delay mismatch for K - 10 and 
for Rayleigh fading, respectively. The ratio of pilot power to data 
power is r - 0.2. We see a difference in the sensitivity to frequency 
errors. Rice fading is more affected by the offset frequency than 

by the Doppler f D . In constrast, Rayleigh fading behaviour is almost 

entirely determined by f D> with little Influence from f Q . 

The irreducible error rate is clearly evident. Its source is the 
phase decorrelation introduced by the delay mismatch during intervals 
of rapid phase change, such as a deep fade. We can obtain a simple 
estimate of the error floor value in the case of Rayleigh fading, the 
most severe condition, by evaluating (17) for small values of f Q r and 

f r . Using the series expansions for sin(), cos(), and J Q (), we 

obtain: 

P eR “ < 2 * f D r)2/8 (18) 

which agrees with Figure 3. To a first approximation, the error floor 
depends only on Doppler, not on frequency offset^ As a numerical 
example, if the error floor is not to exceed 10 , and f^ - 100 Hz, 

then we must keep the delay mismatch |r| < 0.14 ms. Clearly, this is 
not a problem in digital implementations, since this maximum delay is 
only 1/3 bit duration in 2400 bps BPSK. Analog implementations may 
find it a more difficult target. 

A similar approach yields a somewhat poorer approximation to the 
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error floor in the case of Rice fading with values of K over 5 or so. 
Series expansion of the transcendental functions and approximation of 
the integral in the Q(a,b) function yields the error floor: 

P e - K <2„f o ,) 2 e ‘ K (19) 

in the special case of f - 0. 

L 

Narrow PBPF 

One can readily imagine situations in which the combined 
frequency offset and Doppler spread carry part of the received pilot 
spectrum outside the window of the PBPF. Such a situation was 
considered by Larosa & Citron, 1987, for a dual pilot tone system, and 
a number of curves were obtained by numerical integration. However, 
the fading spectrum was taken to be rectangular. Clearly, the 
U-shaped spectrum (2) is more sensitive to band edge distortion 
because a greater fraction of its power is concentrated near ±f D> 

We will model the PBPF as a rectangular LPF with unity gain for 
l f l ^ B p / 2 - With amplitude distortion only, it contrasts with the 

delay mismatch considered earlier, which had phase distortion only. 

The effect on irreducible BER can be determined by evaluating p , 
a, b, and C. We denote the normalized power of the diffuse component 
of the PBPF output by P^: 

p d - Jv f -v | y £ > |2df - Jv f - £ o>v f > df 

- \ arcsin(min(l,(B p -2f o )/2f D )) + ^ (20) 

As for the specular component, we observe that it is present and 
undistorted if it lies in the PBPF passband; otherwise, it is zero. 
Denote by P g the normalized power of the specular component at the 

PBPF output. From (3): 


P 

s 



I Hp (f ) 1 2 df - < 


'1, 

0, 


l f o +f L l * 

otherwise 


Substitution of (20,21) into (5-9) gives: 


B /2 
P 


( 21 ) 


2 

{i ! }-<W(i”A) ! ; C-C l + (22) 

for substitution into the BER expression (26) . 

In the case of Rayleigh fading, we can substitute p from (22) 
into (12), and obtain the error floor: 

P eR “ C 1 • V ~d )/ 2 (23) 

Since should be close to 1 for a reasonable error floor, we can 

apply series expansions for arcsin() and square root to obtain: 

P eR - s C2(£„ * |f c l - B p / 2 )/g 1 '' 2 (24) 
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for small excursions of f Q + | f o | beyond B^/2 . 

As a numerical example, suppose the irreducible BER is to be held 
to 10' 3 . Then from (24), (f D + 1 f Q | B /2)/f D * 8xlO' 5 . This is an 

exceedingly tight constraint; essentially it says that the slightest 
excursion of combined offset and Doppler frequencies beyond the PBPF 
window results in significant performance degradation. 

Bessel PBPF 

Equation (24) suggests that even moderate amounts of rolloff can 
do significant damage if offset and Doppler carry the "horns" of the 
fading spectrum into the rolloff region. Nonrectangular filters, of 
course, do not behave in precisely the same way, and detailed analysis 
requires substitution of their frequency characteristics into (5-7). 

We did this for a sixth order Bessel filter. Of the two integrals in 
(5-7), the one involving |H (f)| was performed analytically; the 

other was evaluated numerically , after a preliminary integration by 
parts to remove the singularity. 

The curves in Figure 3 were prepared with parameters a - f^/R^ 

and P - B /R. . Note that there is an optimum value for the pilot 

filter bandwidth, as we trade noise against distortion of the received 
pilot, and that the optimum depends on K, a, r and E^/N^ The values 

of p were near the optimum. 

For K - 10, we see the onset of error saturation for a - 0.04 and 
P — 0.03; that is, the noise bandwidth of the PBPF is only 0.375 times 
the Doppler bandwidth. However, even for this combination of 
parameters , saturation is a problem only at BERs far lower than we 
normally expect from mobile satellite communication. Rayleigh fading 
is somewhat different. For a = 0.04, we see that the PBPF should be 
at least 1.25 times as wide as the Doppler bandwidth to prevent an 
unacceptable error floor. To put the results in context, oc — 0.04 
represents 96 Hz fading at 2400 bps; that is, it is an extreme case. 
For more reasonable ranges, the sixth order Bessel filter appears to 
be adequate, and it is noise, rather than random FM, which determines 
system performance. 

CONCLUSIONS 

Nonidealites in the pilot tone filter produce an error floor 
which is especially marked in Rayleigh fading with large Dopppler 
frequencies. However, for normal ranges of filter parameters, it is 
not likely to be a problem, and TCT's claim to suppress the error 
floor is justified. 
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Fig. 3. Bessel PBPF, Rice (K-10) and Rayleigh Fading 















